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Abstract 

^ I The photoluminescence in amorphous semiconductors decays according to power law at long 
^ , times. The photoluminescence is controlled by dispersive transport of electrons. The latter is usually 
characterized by the power a of the transient current observed in the time-of- flight experiments. 
Geminate recombination occurs by radiative tunneling which has a distance dependence. In this 
paper, we formulate ways to calculate reaction rates and survival probabilities in the case carriers 
execute dispersive diffusion with long-range reactivity. The method is applied to obtain tunneling 
recombination rates under dispersive diffusion. The theoretical condition of observing the relation 
(5 = a/2 + 1 is obtained and theoretical recombination rates are compared to the kinetics of observed 
photoluminescence decay in the whole time range measured. 
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I. INTRODUCTION 



Photoluminescence in amorphous semiconductors exhibits a power law decay ~ t^^ after 
excitation by weak hght pulses at the absorption edge .^i^i^i'^i^i^i'''i^i^i'^° . The power law decay 
is observed in the long time range from micro seconds to milli seconds. It is widely rec- 
ognized that the photoluminescence originates from the radiative recombination of electrons 
and holes in band-tail states^'^. Holes are quickly trapped and electrons execute hopping- 
random walks among localized band tail states^. At sufficiently high temperatures (T> 
200K), the exponent 6 is close to the value 1.5 predicted from the geminate recombination 
of pairs by normal diffusion However, the exponent is smaller than the value 1.5 at low 
temperatures.— Such deviation is interpreted in terms of deviation of motion of charges 
from the normal diffusio n^i^'^i7i?i?i^? . Murayama and coworkers have systematically investigated 
the exponent 6 by varying temperature. In the experiments, the samples of a-Si:H are excited 
with 10 nsec light pulses of energy around 2.0 eV at the absorption edge.— The photolumines- 
cence main band is observed at the peak energy 1.22 — 1.35 eV and the exponent 6 is independent 
of the excitation intensity in a range from 5[nJ/cm^] to 500[/iJ/cm^].— The measured values 
of exponent gradually decrease by decreasing temperature.— Recently, we have shown that 
the exponent 6 obeys the relation, 6 = a/2 + 1, where the disorder parameter a is the exponent 
characterizing the time evolution of mean square displacements, (r^ (t)) oc a = 1 corre- 
sponds to the normal diffusion and the classical result of 5 = 1.5 is reproduced. The exponent 
a is measured by the time-of-flight technique and the value decreases linearly with decreasing 
temperature . ^i-^^i-^^i-^^i^^i^^i^^ The transient current measured by the time-of-flight technique is 
described by I{t) ~ when t < tr, where a is a constant smaller than 1 and tr is the 

transit time. It is shown based on a continuous time random walk model that a in the transient 
current is the same as the exponent characterizing the dispersive diffusion, (r^ (t)) oc f^M^ 
In the time-of-flight experiments, electrons move into deeper states while they drift from one 
side of the sample to the other. Because electrons in shallow states in the band tails are re- 
leased more rapidly than those in deep states, drift mobility decreases monotonically as time 
passes. The same relation between 6 and a is derived from a similar model using the perco- 
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lation approach to the diffusion of charge carriers in disordered systems.-i^ A simple relation 
between 6 and a, 5 = a/2 + 1, is confirmed by experimental data of a-Si:H taken at various 
temperatures.— The relation between 6 and a expressed as 6 = a/2 + 1 can be obtained from 
a fractional react ion- diffusion equation which is derived from a continuous time random walk 
model.— In the derivation, we have assumed that the recombination takes place at a certain 
distance between the electron and the hole of a pair. However, the most probable mechanism 
of geminate recombination is radiative tunneling. Its rate has a distance dependence, described 
by, koexp (— 2/?r).i^-^ In this paper, we generalize our previous results to include the effect of 
the distance dependence of reactivity. 

Recently, methods to take account of reaction into fractional diffusion equation have been 
developed . '^'''I'^^i^^i^'^i^'^i^^i^^i^'^i^^i^^i^^i^^ Interesting features emerged from such studies are that 
reaction terms have a memory effect because reaction interferes with random walks consisting of 
various hopping frequency associated with distribution in activation energy.— i^SiSi A fractional 
react ion- diffusion equation is derived from a continuous time random walk model by assuming 
that the reaction takes place at a certain distance. In this paper, we formulate ways to calculate 
reaction rates and survival probabilities in the case carries execute dispersive diffusion with 
long-range reactivity. In section 2, we derive methods to calculate reaction rates and survival 
probabilities from a continuous time random walk model with long-range reactivity. In section 
3, we study the square well sink model under dispersive diffusion to confirm the validity of our 
methods by comparing analytical results with those of numerical simulations. In section 4, the 
method is applied to calculate recombination rates in the case of tunneling reactivity of electron 
and hole pairs under dispersive diffusion. Section 5 is devoted to discussion and conclusions. 



II. DISTANCE DEPENDENT REACTION UNDER DISPERSIVE DIFFUSION 



We formulate ways to calculate the survival probability or its decay rate in the case long- 
range reaction proceeds under dispersive transport. We consider geminate recombination of a B 
particle starting at tq, with A. It is straight forward to generalize results for other initial condi- 
tions because the model is linear. B particle migrates by anomalously slow diffusion, oc t" 
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with a < 1. The anomalously slow diffusion is called dispersive diffusion. One of the features 
for dispersive diffusion is absence of characteristic time of jump motion due to trap energy 
distribution. It can be theoretically investigated by a continuous time random walk model in 
which the jumps are performed according to the waiting time distribution . ^^i^^i^^i^^i^^i^^^i^^i^^ The 
waiting time distribution of jump motion in the absence of reaction is denoted by ip (t)- We 
consider ijj [t) which is normalized and has an algebraic asymptotic tail ,^^i^^i^^ 

= ^l^^±iM ^ ^I>±i) , (t>l) (2.1) 

where 7(2, p) = e~*t^~^dt for (Re^; > 0) is the incomplete Gamma function and T{z) is the 
Gamma function.— If we assume activated release, 

^(E) = jrexp(^-^y (2.2) 

and an exponential distribution of activation energy, Eq. 1)2.11) is obtained by 
calculating , 

m= dE-—expi--—)^{E)exp{-^{E)t). (2.3) 
Jo I^bTo \ kBToJ 

a in Eq. ()2.1|) is related to T and To in Eq. ()2.3|) through a = T/Tq. In the absence of reaction 
the above waiting time distribution function leads to dispersive transport, (r^(t)) oc t". On 
the other hand, the waiting time distribution for the reaction in the absence of random walk is 
given by, 

7rc exp [-7,, (r) t] , (2.4) 

where 7rc (r) is the distance dependent reaction rate. In the presence of both jump processes 
and reaction the waiting time distribution of making a jump and that of reaction are given 

jg „10.21. 22.42 

i^outir,t) = V^(t) exp (-7,e (fO ^) • (2.5) 
iprc{r,t) = 7rc (r) exp [-71.C (r) t] / dtitpiti), (2.6) 
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respectively. Eq. (j2.6|) is the reaction rate multiplied by the probability of remaining at the 
site, which decays either by jump motion or reaction. Here, we list the Laplace transforms of 
waiting time distribution functions for later use,^^^*^^^ 



V'out(^s) = ?/'(s + 7rc(r)) 
7rc (r) 



1 - 2F1 [1, a, a + 1, -7r/ s] ~ 1 - 

smTTO; \7r 

s + 7rc (r) \ " 



Tea 



sm Tca 



7r 



V'rc {r, s) 



s + 7rc {r) 



1 - t/) (s + 7rc (r)) 



7rc ir j 



s + 7rc (r) sin vra 



s + 7rc (r) \ ° 



7r 



(2.7) 
(2.8) 
(2.9) 



The equations for the survival probability and the recombination rate can be derived following 
the method presented in a previous paper . ^^1^^ We first formulate the problem on the basis 
of a discrete model on a periodic lattice of dimension d and then take the continuous limit 
in space. We denote the vector characterizing a jump to the nearest neighbor site j, by bj 
(j = 1, 2, ■ ■ ■ , 2d) and the jump length by b. The equation for the probability ri{fi, t) dt of just 
arriving at site rj in the time interval between t and t + t is written as, 

'7(^*5^) = 771 / dtiilJout{ri-bj,t-ti)ri{fi-bj,ti) + 6f^^ro- (2-10) 
^" j^i Jo 

By subtracting Jq dti'tpoutiri,t — ti)r]{fi,ti) from both sides of Eq. ()2.1()|1 . we obtain in the 
small b limit, 

b^ 



1 - V^out(^, s) 17 (r, s) = y V^V'out(^, 5)57 (f, s) + 5{f-fo) 



(2.11) 



where the Laplace transform is introduced, i.e., fj{f, s) = dtexp{—st)ri{r, t). Eq. ()2.1H] can 
be rewritten as, 

62 



1 - i){s) fi{f, s) = —V ipout{r, s)fi (f, s) - 6ip{f, s)fi (f, s) + 6{f- fo) 
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where 6ijj{f, s) = ip (s) — '?/'out(^, s). Up to this point b has been assumed to be small but finite. 
In the continuous limit, we have to take the limit of 6 — ^ with the generalized diffusion 
constant given by^ 



smvro; 
27ra 
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au2 



(2.13) 



kept constant. Therefore, we should consider the hmit, 7rc(r)/7r ^ which leads to 
V'out {r, s) — > 1 for s ^ 0. Since iprc {r, 0) = 1 — -ipout {r, 0), which follows from the fact that par- 
ticles at a given site perform either jump or reaction, i.e., V'out {r, t) dt + iprc {r,t) dt = 1, 
we find, 

6^^(f,s)^^Prcir,0)^^^(-^^^'] fors^O. (2.14) 

smvra V 7r / 

With the aid of the above relation together with Eq. ()2.7|) . Eq. ()2.12j) becomes in the continuous 
limit as, 

s"fi{f, s) = DaV'^fj (f, s) — ka (r) fj (r, s) + ^^^^"^ 7" (5 (f — tq) , (2-15) 
where the generalized reactivity is defined by. 

The effective reactivity ka (r) depends on a which characterizes the dispersive diffusion. The 
dispersive diffusion competes with the reaction kinetics. As a result, the characteristic exponent 
a enters into the expression of the effective reactivity. The recombination rate R{t) satisfies, 

R{t) = dhj dfij,,{f,t-h)7]{f,ti). (2.17) 

After the Laplace transform, it is expressed in the asymptotic limit of s ^ as, 

R{s) = j dfip,c{f,0)ri{f,s) . (2.18) 

The recombination rate R{t) and the survival probability N{t) are given by, 

R{t) = -§-^Nit) = j df'ip,,{f,0)r]{f,t), (2.19) 

N{t) = j df^,,{f,0) dh7]{f,t,). (2.20) 

Therefore, survival probabilities and recombination rates are expressed in terms of the prob- 
ability T] (r*, t) dt of just arriving at site f in the time interval between t and t + dt, which is 
obtained by solution of Eq. ()2.15p followed by the inverse Laplace transform. 
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Although recombination rates and survival probabilities can be calculated by the above 
method, there is more convenient alternative method which naturally reduces to the con- 
ventional method in the case of normal diffusion, a = 1. By introducing, p{f,s) = 
1 — ipi^s) 1] (r, s) /s, Eq. ()2.15|) is rewritten as. 



sp (f, s)-5{r- fo) = " [DaV^p (f, s) - K {r) p (r, s)] . 
After inverse Laplace transform, a fractional reaction-diffusion equation is obtained. 



(2.21) 



d_ 
di 



p (f, t) 



d_ 

di 



dti 



r(«) it - h 



— [D^V^p (f, h) - K (r) p (r, ti)] . 



(2.22) 



Eq. (|2.18|) is transformed into. 



R{s) = s 



l-a 



dfka (r) p (r, s) 



(2.23) 



by noticing the relation, p{f,s) = 1 — ip{s) r]{f,s) /s. After the inverse Laplace transform, 
Eq. (IT^ becomes. 



R{t) = -|iV(t) 



d 



dti- 



vl-a 



ka (r) p (r , ti 



(2.24) 



dtj, ■^"'T{a){t-h) 

In principle, the density pd (^5 1) should be given in terms of the probability of remaining 
at a site, (f){r,t) = /^°° rfti [^'rc (r, ti) + V'out (r, ti)], as Pd{r,t) = J^^ dticp {r,t - ti) rj {f,ti). 



After the Laplace transform, we obtain pd {r, s) 



1 - i'out {r, s) - iJrc (r, s) r]{r,s) /s. 



The Laplace transform of the density is expressed in terms of p (r, s) as pd {r, s) = 
1 ~ 4'ont (^) s) — iprc s) p{rjS)/ 1 — ip{s) . In the case of normal diffusion, a = 1 holds 
and the above expression for the density reduces to p (r, s) as shown in Appendix El but they 
are not equal for a < 1. Eq. ()2.24j) is derived from Eq. ()2.18|1 . On the other hand, if we 
assume N{t) = J dfp{f,t), integration of Eq. ()2.22|) over r also yields Eq. ()2.24j) . There- 
fore, the relation N{t) = J dfp (r, t) should hold, although p (r, t) without integration over 
r cannot be interpreted as the real density for a < 1. Eq. ()2.24|) together with Eq. ()2.22p 
is an alternative method to calculate recombination rates and survival probabilities. They 
are useful because recombination rates and survival probabilities under dispersive diffusion are 



derived from those under normal diffusion with a = 1 by substitution, Di 



Dn and 



ki (r) s^^°'ka (r) after Laplace transform. In the following section, we study the validity 
of such approach by comparing with numerical simulations. It should be also noted that the 
effective reaction rates depend on a through Eq. ()2.16p . In the case of tunneling with rate given 
by 7rc [r) = koe'"^^^, the effective reactivity ka (r) is given by ka (r) = /cqC"^"^''. The effective 
tunneling distance 1/ (a/?) increases when a is decreased by energetic disorder. In the case of 
energy transfer, 7rc (^) = ko^Rp/r)^, with Forster radius Rp, the effective reactivity is given 
by ka (r) = k^lRp/r)^"'. The exponent 6a is decreased with decreasing a and the long-range 
nature of energy transfer is enhanced. Shushin obtained similar results. He showed that the 
effective reaction radius appearing in the steady state solutions of Smoluchowski-type stochas- 
tic Liouville equations increases with decreasing a for tunneling reaction.^^' He also showed 
that the kinetics changes at a = 1/2 for energy transfer of Forster type.— We do not further 
investigate the energy transfer but focus our attention on tunneling recombination reactions. 



III. SQUARE WELL SINK MODEL UNDER DISPERSIVE DIFFUSION 



In this section, we consider the model in which reactivity constant is k^ within the sphere of 
radius R.4S4£iiLi£ In the absence of diffusion, the waiting time distribution of reaction is given 
by Eq. (jOI) with 

7re(r) = A;o*S(r), (3.1) 
where the reaction function S{r) is given by. 



S 



r) 



1 ioT r < R 
for r > R. 



Reaction takes place with the uniform reactivity whenever the distance between reactants is less 
than R. This model was already studied in the case of the normal diffusion with a = i.i^iiSiiLI^ 
The effective reactivity is obtained from Eq. ()2.1(ij) as 

ka{r) = k^*S{r). (3.2) 
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The effective reactivity is scaled by the exponent a < 1 under the dispersive transport. As 
described in the last part of the previous section, we apply substitution, Di s^^°'Da and 
ki (r) — s> s^~'^ka {t) to the results derived for normal diffusion in the Laplace domain. We 
consider the geminate recombination in which the initial distance between the pair is tq. The 
escape probability, ip {tq) = N (oo) is obtained by substituting ko/Di — > k^/Da into Eq. (5.13) 
of Tachiya's result,— 

D / tanh (R,/k^/D^)\ 
y,(ro) = 1-- 1 \ M . (3.3) 

The Laplace transform of the survival probability is obtained from Eq. (5.10) of Tachiya's 
result,— 



1-1-^ 



RJ'H^R cosh fl, M - s,uh rJ''° + *» 



(3.4) 



The numerical Laplace inversion of Eq. ()3.4|1 by Stehfest algorithm^ is employed to get the 
theoretical results. We compare these results to those of numerical simulations. 

The simulation procedure is the same as that used previously.— We explain it briefly here. 
In our model one reactant is fixed at the origin, and the other performs a random walk starting 
from an initial separation, ro = 2R. The random walk is realized as a sequence of detrapping and 
instantaneous trapping events, with the trap energies generated according to the exponential 
distribution characterized by an attenuation parameter /cb^o- The rate of release from the 
traps is assumed to be given by Eq. ()2.2j) . and the detrapping time from a trap of energy 
E is obtained from the exponential distribution with the mean value l/'-f{E) . The Gaussian 
distribution of the jump length is assumed for each direction of the Cartesian coordinates with 
the standard deviation of b and the zero mean value. We assume that the reactions occur within 
a spherical shell of radius R and the reactivity is given by /cq- Iii the simulation, the reactions 
are modeled in the following way. When a reactant is in the reactive sphere, we generate the 
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reaction time tj-c from the exponential distribution with the mean value 1/ fco ? and compare it 
with the detrapping time t calculated for the current trap. If trc < t then the reaction occurs, 
otherwise the particle jumps to another trap and the simulation continues. The trajectory is 
calculated until the particle either reacts or escapes to a large distance rmax = lOOR . The 
simulation is repeated for at least 10'* independent trajectories. 

The calculated diffusion coefficient is in good agreement with Eq. (j2.13|) . The simulation 
results of the escape probability against the reactivity are shown in Fig. 1. Excellent agreement 
is found between the analytical result of Eq. (jH.Hj) and simulations regardless of a values. In 
Fig. 2, the survival probabilities obtained from the inverse Laplace transform of Eq. ()3.4p are 
compared with the results of numerical simulations for a = 0.5. Theoretical results coincide 
with the simulation data for the chosen values of reactivity. These results confirm the validity 
of procedure to calculate the survival probability and escape probability in sub-diffusive media 
with long-range reactivity described in the previous section. In the next section, this procedure 
is applied to the case of the tunneling recombination under dispersive diffusion. 



IV. TUNNELING RECOMBINATION UNDER DISPERSIVE DIFFUSION 



In this section, we examine recombination rates under the radiative tunneling described by 
the first order rate, koe~^^'^, where (3 is known to be ~ 0.05 — O.lA ^.^ The effective reactivity 
under the dispersive transport is given by, ka {r) = exp (— 2a/5r), and we have to solve, 

sp (r, s) - ^^p^^ = s'~^ [D^V'p (r, s) - k^ (r) p (r, s)] . (4.1) 

The kinetics of geminate recombination by normal diffusion and tunneling was investigated 
previously. However, a perfectly absorbing boundary condition was imposed on a sphere of 
atomic dimensions at the origin with a certain radius smaller than the initial distance of a 
pair of an electron and a hole, In this paper, we assume recombination occurs only by 

tunneling and impose a perfectly reflecting boundary condition at the radius, R. By introducing 



transformation, p (r, s) = a {r, s) /r , and a new variable, y = ^/k^^JD^exp {—a/3r) / {a/3), the 
homogeneous part of Eq. 1)4.11) becomes the modified Bessel differential equation . 
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Thus, the homogeneous solutions of Eq. (j4.H) are, 



P{r,s) = ^[l^J^e-^^], (4.2) 



4{r,s) = -K, J-./^e--/^M. (4.3) 

r yap \ Da J 

General solutions are written by introducing constants, Ci, C2 and C3, which are determined 
to satisfy the boundary conditions. 



P [r, s) 



Cip (r, s) r > ro 

C2q (r, s) + Csp (r, s) r < tq 



Solution should be continuous at r = ro and satisfy Eq. ()4.1|1 . 

CiP (ro, s) - C2q (ro, s) - C3P (ro, s) = (4.4) 
C2^J (ro, s) + (ro, .) - C,^P (ro, .) = (4.5) 

Perfectly reflecting boundary condition at r = < ro leads to, 

C2Q^q {R, s) + C,—P {R, s) = 0. (4.6) 

Later we take the limit of R —>■ 0. We have 3 equations for 3 unknown coefficients, Ci, C2 and 
C3. By solving Eqs. ()4.4|) - ()4.6|) . the solution in the Laplace domain is found, 

-r-F] — n Q {ro, s) - p ro, s) ]p{r,s) r> ro, 

'g(r,g)- j p(r, g) hp(ro, g) r < ro. 



47rL)„a/5 ' ' ' P'{,R,s) 

where pf{r,s) = -^p{r, s) and qf{r,s) = -^q{r,s). In the above derivation, the relation, 
p{r, s)q f{r, s) — p f{r, s)q{r, s) = aP/r"^, is introduced. Because the sample is excited with 
relatively weak light pulses at the absorption edge, an electron is created in the vicinity of a 
hole. The solution is simplified when ro = R, 

AnDaR pf [R, s) 
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which is further simphfied in the hmit of i? — > as, 



-1 ^A,.l£L al 



P [r, s) 



al3y Da 



aP V Da 



1 / sQ 



al3y Da ^ 



(4.9) 



Laplace transform of the survival probability is calculated by, N (s) = An r p (r, s) d r, and 

Jo 

the result is expressed as. 



N(s) 



«-i \ 2ap V D. 



Uc 



1 / 

al3\J Da 



E 



2a/3 V Da, 



2n 



Da \ a(3 V Da J 



fiV Da) {"'^ afiV Da 



(4.10) 



where the series expansion of the modified Bessel function is introduced,— 



2n 



n=0 



n\T{u + n + 1)' 



By taking the limit of s 0, the escape probability is obtained as, 

1 



ip{ro = R-^0) 



1 



a/3 V D. 



(4.11) 



(4.12) 



In the case the intrinsic recombination rate is small, tt^a/ tt- < 1, the recombination rate in 

' zap y Da — ' 



the long time limit becomes, 

dN{t) 



m 



akn 



1 



dt 



r(l-|) {2a(3^/K) (l + 



ftp 



2 * 



(4.13) 



as described in Appendix |B] Clearly, the relation, 5 = 1 + f , holds. 

In the limit of 2^\/ if" recombination rate in the long time limit is simplified as. 



R{t) 



akn 



r(l-f) {2a(3^y 



ta/2+l ' 



(4.14) 
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Therefore, the recombination rate depends on the disorder parameter a in the form, 

m - (4.15) 

where C = 7r'^^/fco is a constant. The relation, 5 = 1 + ^, is observable when is small 
compared to the rate of escape of an electron from the region within the distance ~ 
10 — 30A from a hole by hopping. In this case an electron may enter that region many times 
before recombination. In the previous paper, the measured photoluminescence is shown to obey 
a scaling relation, R (t) ~ ^q^q/2+i ' "where ~ 6.8 x 10^^[s~^].^'' This value of is obtained 
from the temperature dependence of photoluminescence intensity observed at a delay time of 
10[//s].— The temperature dependence is assumed to arise from the disorder parameter, which 
is phenomenologically described by a linear function, 

T 

a = ao + —, (4.16) 

where ao = 0.15 and To = 313.4 [K].i2*i2*i^ If we assume an exponential distribution of band 
tail states characterized by an attenuation parameter /cb^o? is given by a = T/Tq. The 
presence of the constant term, oq, in Eq. ()4.16p implies that activated release of an electron 
from an exponentially distributed band tail states is an oversimplified model. Nevertheless, the 
linear temperature dependence is a good approximation for the observed values of the disorder 
parameter, a. Eq. (j4.14p also has the scaling form with a constant, C = 7r^^/fco- Although 
radiative tunneling rate has a distance dependence, photoluminescence can be modeled by 
recombination reaction at encounter distance as we have done in the previous paper, as far 
as long time asymptote is concerned. The tunneling rate competes with the jump frequency 
of electrons only at localized sites in the vicinity of a hole. The specific condition is given 
by 2^\J ^ < 1- Note also that the inverse /? of tunneling distance in the expression of 
photoluminescence decay scales by a, i.e. , a(3. In other words, the effective tunneling distance 
increases when a is decreased by energetic disorder . In Fig. 3, theoretical lines obtained from 
the inverse Laplace transform of -R(s) = 1 — siV(s), where N [s) is given by Eq. ()4.10|) . are 
presented together with the measured photoluminescence decay at 145 [K] 71. is varied while 
keeping i^r^"^ /k^ ~ 6.8 x 10^^ [s~^]. This value of ^7r^^/fcoj guarantees the amplitude of 
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asymptotic power law decay at other temperatures as explained above. The long time asymptote 
is well approximated by power law decay, t^^ with 5 = 2/a + 1, where a ~ 0.6 at 145 [K] is 
obtained from Eq. ()4.16p . The asymptotic time region described by power law increases with 
increasing 7r. In the case of 7r = 10^^ [s~^], the line in Fig. 3 is described by power law decay in 
the entire time region of the Figure. Experimental data are close to the line of 7r = 10^ [s^^]. In 
the case of 7i. = 10^ [s~^], the onset time of asymptotic power law decay appears later than that 
of experimental data. Before the onset time, the slope increases with increasing the delay time 
after the light pulse. The value of 7r = 10* [s~^] is smaller than the values 7r ~ 10^^ — 10^^ [s^^] 
estimated from the time-of-flight technique.— i^^*^ In Fig. 4 7i. is fixed to the value of 10^^ [s^^] 
and the value of ko is varied. The relation ^7r''^/fcoj ~ 6.8 x 10^^ [s~^] holds for ~ 10^^ [s^^]. 
For other values of fco, the temperature dependence of photoluminescence intensity at 10[yus] 
is not reproduced. Nevertheless, we look for the possibility of fitting experimental values of 
photo-luminescence at 145 [K], by assuming 71. = 10^^ [s^^]. In the case ko is smaller than 71., 
namely, ^tpV Id' — photoluminescence decay is described by power law, with the 

exponent 6 = a/2 + 1 ~ 1.3 in the entire time region of the Figure. In the case of ko > 7r, 
namely, -^^^J^- > 1? the photoluminescence decay is described by power law with the 
exponent larger than the value, a/2 + 1 ~ 1.3. Since the experimental data deviate from the 
power law early times, they are not reproduced at early times by theoretical lines 

with 7i. = 10^^[s^^]. Only the asymptotic decay is reproduced in the case of -^sj^ < 1- 
When fco is large, namely, 2^\j if" > 1; the photoluminescence decay is not described by the 
algebraic decay expressed by Eq. ()4.14j) in the time range described in Fig. 4. Nevertheless, 
the decay is approximately described by the power law of time with the exponent larger than 
a/2 + 1. In this limit, not only hopping but also the distance dependence of the intrinsic rate 
controls the main part of the photoluminescence decay. 



V. DISCUSSION AND CONCLUSIONS 



We have derived methods to calculate reaction rates and survival probabilities for a con- 
tinuous time random walk model with long-range reactivity. In Eq. ()2.19|l - ()2.20j) . survival 
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probabilities and recombination rates are expressed in terms of the probability of just arriving 
at r in the time interval between t and t + dt, which is obtained by the inverse Laplace trans- 
form of the solution of Eq. ()2.15p . We also present more convenient alternative method, Eq. 
fl2.24|) together with Eq. ()2.22|) . In this method, recombination rates and survival probabilities 
under dispersive diffusion are derived from those under normal diffusion with a = 1, by substi- 
tuting Di —>■ s^^°'Da and ki (r) s^~°'ka (r) in the Laplace domain. The effective reactivity is 
related to the intrinsic reaction rate by Eq. (|2.16p . The effective reactivity in the case of tun- 
neling, 7rc if) = /cqc"^^'', under the dispersive diffusion is expressed as ka (r) = /cq exp (— 2a/5r). 
Therefore, the effective tunneling distance increases when a is decreased by energetic 

disorder. In the case of energy transfer, 7^ (r) = ko^Rp/r)^, with Forster radius Rp, the effec- 
tive reactivity is given by ka (r) = fcg (i^i?/''")^"- The exponent 6a is decreased with decreasing 
a and the long-range nature of energy transfer is enhanced. 

In order to confirm the validity of the methods, we study the square well sink model un- 
der dispersive diffusion . ^^'^i^^'^'i''^^ The escape probability and the Laplace transform of survival 
probability are known for normal diffusion.— Therefore, these quantities under dispersive dif- 
fusion are obtained by introducing substitution in the Laplace domain, Di s^~°'Da and 
fco — > s^~"A;q . The results are compared to direct numerical simulations of random walk and 
excellent agreement is found between them. 

Finally, the method is applied to obtain recombination rates by tunneling under dispersive 
diffusion. In the previous work, the relation 5 = a/2 + lis confirmed for the experimental data 
of a-Si:H and the scaling relation of the amplitude, R{t) ~ ^c»^i/2+i with ~ 6.8 x 10^^ [s~^] is 
obtained.— We found that these relations are observable in the case k^ is small compared to the 
rate of escape of an electron from the region within the distance 1/ {a(3) ~ 10 — 30A from a hole. 



namely, 2^ y 'd' "^^^ experimental data is close to the theoretical line with 7r = 10^ [s^ ] 
over the whole time range measured. However, the value of 7r = 10^ [s~^] is smaller than the 
values 7i. ~ 10^^ — 10^^ [s~^] estimated from the time-of-flight technique.— ^^^^ The discrepancy 
could be due to some unknown factors which may change the recombination rate such as 
the initial distribution of an electron around a hole, electrostatic interaction energy between 
them. In the presence of distributed site energies of localized states, Coulombic interaction 




15 



can be ignored in the case the distance between a hole and an electron is larger than the 
characteristic length Rs = ^/e%/ (eA;B7o) ~ 40A, where e ~ 10 is the dielectric constant of 
a-Si:H,^° This distance is relatively small but still it is comparable to the hopping distance. It 
should be also noticed that the hopping frequency estimated from the time-of-flight technique 
is also approximate. Since the correct value of the hopping frequency is not known, it is not 
meaningful to introduce the distribution of initial positions of electrons which is not known 
either. When ko is large, namely, 'd~ ^ photoluminescence decay is approximately 

described by the power law of time with the exponent larger than a/2 + 1 as shown in Fig. 4. 
In this limit, not only hopping but also the distance dependence of the intrinsic rate controls 
the main part of the photoluminescence decay. 
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APPENDIX A: RELATION BETWEEN p (r, s) AND THE LAPLACE TRANSFORM 
OF THE REAL DENSITY 



By using the first equalities of Eqs. ()2.8p and ()2.9p . we have. 



1 - V'out {r, s) - V'rc {r, s) 



1 - (s + 7rc (r)) 



s + 7rc (r) L 

Substituting the second equality of Eq. ()2.8|) into the above equation yields 



(Al) 



1 - V'out {r, s) - tprc (r, s) 



s + 7rc (r) 



(A2) 



s s + 7rc [r) sin na \ 7r 

In the case of a = 1, the right hand side of Eq. ()A2|) is independent of 7i.c (r) and can 
be shown to be identical to 

fj (r, s 



r]{r,s} 



l—4'{s) /s by Eq. ()2.7|) . Therefore, p{r,s) defined by 
1 — ■?/' (s) /s coincides with the Laplace transformation of the real density given by 
1 ~ ^out — ipj-c s) / s. However, they are generally different for a < 1. 
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APPENDIX B: DERIVATION OF EQ. 



By substituting Eq. (j4.11|) into Eq. (j4.1Uj) . we obtain, 



N{s) 



V a/ 



2n 



In the limit of s ^ 0, Eq. ()B1|) leads to 



(Bl) 



R{s) = l- sN{s) 



2aP V D. 



1 k, 



aP V D, 



(B2) 



The inverse Laplace transform using Tauberian theorem yields Eq. (|4.13|) . By expansion of 
Eq. ()B1|) in terms of series of and then performing the inverse Laplace transform, we 

obtain the survival probability. 



N{t) ~ EiL 



aP\/ Da + 



^0 ] tc./2 



+1 



(B3) 



expressed in terms of the generalized Mittag-LefHer function defined by^ 



k=0 



r {ak + h) ' 



where Ea {z) = Ea^i {z). In the same way the recombination rate is obtained as 



R{t) ~ 1 - ^£;|,o 



a(3\/Da + 



t 2'2 



iapy/Da 
aP\/Da + 



.a/2 



ua 



.a/2 



(B4) 



(B5) 
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We can confirm that the above equation reduces to Eq. (j4.13|) by the asymptotic expansion of 
the generalized Mittag-Leffler function,— 

-k 



z 



In the case of normal diffusion, the inverse Laplace transform of Eq. ()B1|) becomes 

2 



(B6) 



N{t) 



1 



1 + 



1 , ^0 
1 + . ^ exp 



4/32 Di 



t 



erfc 



/3VDi + 



(B7) 



41 



where erfc(2;) = ^ dtexp {—t^); 

In this paper, we focused our attention on the limit of tq = R. As a reference, we list a 
result for tq > R. In the case of i? 0, Eq. ()4.7j) is rewritten for an arbitrary initial position 

vq as. 



P[r>,r<,s) 



where r> = max (r, ro) and r< = min (r, tq). 



g |^r<, s) 7 , V P [r<, s) 



\ 



I. -I — 

C I a/3 V -Da 



p(r>,s), (B8) 
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Fig. 1: The escape probability versus 0^/ {R^kQ ). The initial distance is tq = 2R and 
b/R = 0.1. The symbols denote simulation data. The circles, triangles and crosses correspond 
to the a values of 0.5, 0.75 and 1.0, respectively. The line shows the analytical result of Eq. 

Fig. 2: The survival probability as a function of time for a = 0.5. t is normalized by 7r, 
ro/R = 2 and b/R = 0.1. ko {R^/D^f'"" = 10°, 10\ 10^ from top to bottom. Sohd lines are the 
numerical Laplace inversion of Eq. (j3.4j) . Dots represent the results of numerical simulations. 

Fig. 3: Decay of photoluminescence in a-Si:H excited with the pulsed laser of photon energy of 
2.32 [eV]. Dots show the observed decay in the whole time range measured. The photolumines- 
cence has been observed at 145[K].iS*iS Lines are calculated from the inverse Laplace transform 
of ^(s) = 1 - sN{s), where N (s) is given by Eq. (ICTID with (^Tr^V^^o)^ ~ 6.8 x lO^^is"^] 
being fixed, a = 0.15 + 145/313.4 ~ 0.61, b = lOA and /3 = 0.05[1/A]. Amplitudes are adjusted 
by comparison to experimental data. The short dashed line, solid line and long dashed line 
represent 7r = 10^ [s~^], 7r = 10^ [s""*^] and 7r = 10^° [s~^], respectively. 

Fig. 4: Decay of photoluminescence in a-Si:H excited with the pulsed laser of photon energy 
of 2.32 [eV]. Dots show the observed decay in the whole time range measured. The photo- 
luminescence has been observed at 145 [Kj.^^^ Lines are calculated from the inverse Laplace 
transform of -R (s) = 1 — sN{s), where N (s) is given by Eq. ()4.10p with 71. = 10^^[s~"'^], 
a = 0.15 + 145/313.4 ~ 0.61, b = lOA and P = 0.05[1/A]. Amplitudes are adjusted by compar- 
ison to experimental data. The thin solid line, thick solid line and short dashed line represent 
ko = 10® [s~^], ko = 10^"*^ [s and ko = 10^^ [s "*^], respectively. They obey almost the same power 
law, l/t^, with the exponent 6 = a/2 + 1 ~ 1.3. The dashed-and-dotted line and long dashed 
line represent ko = 10^'^ [s~^] and ko = 10^^ [s^^], respectively. They are well approximated by 
the power law, ~ 1/t^'®, with the exponent 1.6 larger than 6 = a/2 + 1 ~ 1.3. 
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Fig. 1, K. Seki, M. Wojcik, and M. Tachiya 
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Fig. 2, K. Seki, M. Wojcik, and M. Tachiya 
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Fig. 3, K. Scki, M. Wojcik, and M. Tachiya 
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Fig. 4, K. Seki, M. Wojcik, and M. Tachiya 
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